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We propose an optical lattice scheme which would permit the experimental observation of Zitterbewegung 
(ZB) with ultracold, neutral atoms. A four-level "tripod" variant of the usual setup for stimulated Raman adi- 
abatic passage (STIRAP) has been proposed for generating non-Abelian gauge fields Dirac-like Hamil- 
tonians, which exhibit ZB, are simple examples of such non-Abelian gauge fields; we show how a variety of 
them can arise, and how ZB can be observed, in a tripod system. We predict that the ZB should occur at 
experimentally accessible frequencies and amplitudes. 



A driving force behind the study of ultracold atoms is their 
potential use as highly tunable quantum simulators for physi- 
cal systems, ranging from quantum phase transitions in solids 
1^21 to black holes | 3 1. In particular, the high degree of control 
over length and time scales in cold atom experiments allows 
for the possibility of observing phenomena that are experi- 
mentally inaccessible in their original counterpart systems. In 
this paper, we propose an experiment which simulates the rel- 
ativistic (and recently, controversial |4 |) phenomenon of zit- 
terbewegung (ZB), a jittering motion caused by interference 
between the positive and negative energy components of the 
wavef unction of a Dirac fermion. 

For a relativistic electron, the ZB frequency is of the order 
of mc^/h ^ 10^^ s~^ and the amplitude comparable to the 
Compton wavelength, h/mc ^ 10~^^ m. ZB has never been 
observed for free electrons, as these time and length scales 
render it experimentally inaccessible. The presence of ZB is, 
however, a general feature of spinor systems with linear dis- 
persion relations. Trapped ions Q as well as condensed mat- 
ter systems, including graphene (61 |7l El and semiconductor 
quantum wires ||9j[T0l, have been proposed as candidate sys- 
tems for observing ZB. 

In this paper, we propose a scheme for observing ZB in 
ultracold neutral atoms. A four-level "tripod" variant of the 
usual setup for stimulated Raman adiabatic passage (STIRAP) 
has previously been proposed for generating non-Abelian 
gauge fields 1 1 1. Dirac-like Hamiltonians, which exhibit ZB, 
are simple examples of such non-Abelian gauge fields, and we 
show how a variety of them can arise in a tripod system. The 
Hamiltonian for atoms in an optical lattice is Dirac-like in the 
subspace of the tripod's two degenerate dark states. We pre- 
dict that an atom's mean position should thus undergo Dirac- 
like ZB. However, the characteristic amplitude of tripod ZB 
is the optical lattice wavelength, vs. the Compton wavelength 
of Dirac ZB, and the oscillation energy is proportional to the 
lattice recoil energy vs. the rest mass energy of the electron. 
This places tripod ZB well within the range of experimental 
observation, with a characteristic frequency of MHz vs. the 
THz domain predicted for condensed matter implementations 
ISl . Although we treat here the case of a noninteracting gas, 
the ZB persists under the addition of weak interactions: The 
Hamiltonian separates into center of mass and relative coordi- 
nates, and the center of mass Hamiltonian is again Dirac-like. 
In a dilute atomic cloud, the ZB should thus manifest itself as 



an oscillation of the cloud's center of mass. 

We consider the tripod STIRAP scheme described in |[T| 
and shown in Fig. [T] The Hamiltonian in the interac- 
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Figure 1: Left: Tripod STIRAP scheme with Rabi frequencies as 
defined in Eqs. Right: Schematic geometry of tripod laser 

beams yielding the Rabi frequencies of Eqs. 

tion picture is H = —hY4=i^i\0){i\ + h.c. Defining 0/ = 
\jY!j=\ l^il^' the dressed states include two dark states de- 
generate at zero energy: \D\) = ^ (^l2|l) — |2)) and 

\D2) = ^ (^1J^13|1)-^^2^3|2)-0^|3)) (we have chosen 
an orthonormal basis). Suppose the atoms are now slowly 
moving in the field. The degeneracy causes the Born- 
Oppenheimer approximation to break, yielding an effective 
U(2) non-Abelian gauge field. The effective Hamiltonian in 
the 2 X 2 dark subspace is H = ^{p — A)^ -\- ^ where m is 

the atom's mass, Aij = ih{Di\V\Dj) is an effective vector po- 
tential, and <i> is a scalar Born-Huang potential resulting from 
the coupling to the bright subspace. The following choice of 
Rabi frequencies 
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for Q,, 8, and ko constant corresponds to the laser beam ge- 
ometry in Fig. [TJb), and yields a Dirac-like Hamiltonian (a 
related setup was proposed in 1 11 1 in the context of observing 
spin relaxation effects). Specifically, after some trivial gauge 
transformations, the two dark states feel an effective vector 
potential 
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and an effective scalar potential <E> = VoG^^ 
h kr, /, _9\ X " [Ay.Az] 0, and the field is non- 



where Vq = 

^ (l -£^) . In general, 

Abelian. The vector potential in Eqs. (|4][5]) is valid for 
— 1 < 8 < 1, acquiring an extra normalization otherwise. 
We can write the full Hamiltonian in the dark subpace as 



H = 
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[(1 - e^)pya, - 2Ep,ay] -^VqG,; (6) 



we shall henceforth not consider the free x direction. Think- 
ing of |/)i(2)) as "spin up (down)," the Hamiltonian of Eq. ^ 
is effectively a spin-orbit coupling, in the presence of a homo- 
geneous magnetic field along the z direction. In fact, letting 
a = ^8, defining 8/? = — 1 + a/2, and choosing 8 = Er, we 
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have H = ^ -\- a{pyGz — Pz<^y) -\- VoG^-thc Hamiltonian for a 
2D electron gas in the y-z plane, with Rashba spin-orbit cou- 
pling, and subject to a homogeneous magnetic field along the 
z axis. Ref. 1 12] proposes an alternate scheme for generating 
spin-orbit coupling with ultracold atoms. 

It is possible to remove the scalar potential by applying a 
state-dependent external potential to the system. Denoting the 
potential felt by |/) as Vi{f), choosing Vi(r) = V2{r) =V{f), 
andy3(r) =y(r)+Vo(l+8^)/(l-8^) subjects the dark states 
to an additional potential V = V (r) (g) I — <!>. If the scalar po- 
tential is thus removed from Eq. ([6]), the resulting Hamilto- 
nian is extremely versatile, for two reasons: (1) 8 is tunable 
and (2) the dark states form a degenerate subspace, for which 
any basis is equivalent. In fact, depending on the direction we 
call "spin up," this Hamiltonian can be viewed as a variety of 
Dirac-like Hamiltonians (see Table |l|. 

The eigenstates of the system are spinors of the form 
e^^y^e^^^^^li'Ji) with / = a^b. The dispersion relation of the 
Hamiltonian in Eq. ([6]) is 




Figure 2: (Color online) Energy surfaces in momentum space for 
^^Rb with ^0 = 27i/820nm~^ and 8 = 8/?. A conical intersection oc- 



curs at {ky.i^z) = V^o 



kz) = (^0'^)' ^ circuit of this conical intersection in k 
space gives a Berry phase. 



expansion coefficients to be k dependent (and also multival- 
ued). The presence of the Bem^ phase thus translates into a 
k dependence of ^(^) in Eq. (jsj); we later show that it is this 
nonvanishing V^^(^) that gives rise, directly, to ZB. 
Considering an initial wavepacket 
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it is straightforward to show that its time evolution is 
\|/(r;f) = ■^jdkg(k-())e"''e-''^^^^^'l'' 
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The energy surfaces in Eq. ^ have a conical intersection at 
2]). A circuit of the degeneracy in 
erry phase: Defining 



{k^X) = {k^.O) (see Fig. 
momentum space yields a 
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\b;k) 



iCOS 



-sin^y^l^ which 



are multivalued for a particular (ky^k^). The momentum- 
space Berry phase is much like the one encountered in 
graphene, which gives rise to phenomena like the half-integer 
quantum Hall effect [14J. We now examine the role that the 
Berry phase plays in generating ZB. 

Consider the time evolution of a Gaussian wavepacket pre- 
pared in a superposition of dark states. In the non- Abelian 
case (0 < 8 < 1), the eigenvectors have an associated Berry 
phase, and are both k dependent and multi-valued. How- 
ever, the initial spin state must be single valued, forcing its 
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where co(^) = ^{E^{k) -E-{k)). 

ZB, in the Dirac equation, is an oscillation of the average 
position (r(^)). The usual method of understanding the Dirac 
equation, and related equations is to derive equations of mo- 
tion for the Heisenberg operators, and show that they oscillate 
in time |71[8)|9)110|. We instead work in the Schrodinger pic- 
ture, which makes explicit the connection to Berry phase. As 
[H^p] = 0, it is convenient to work in the momentum basis, 
and calculate 



m)=ijc 



dk^\k\t) -VMku) 



where \j/(r;r) = / dk^(k\t)e^^'^ / {2%) is the spinor wavefunc- 
tion. After some algebra, we find 
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Graphene sheet, in vicinity of Dirac point 



Table I: By tuning 8 and choosing different states to represent spin-up, the tripod setup can replicate a variety of Dirac-like Hamiltonians (see 
e.g. L2J). Ref. I.13J proposes an alternate method of generating a graphene-like Hamiltonian. 



where the final term, which oscillates in time, is ZB. The am- 
plitude of the oscillation is proportional to V^^(^). We had 
previously shown that the k dependence of i^Qc) occurs as 
a direct consequence of the eigenfunctions being multival- 
ued. The Schrodinger picture thus illuminates what is not evi- 
dent in the Heisenberg representation-that the ZB here can be 
viewed as a measurable consequence of the momentum-space 
Berry phase. 

We now suggest a possible experimental demonstration of 
ZB using ultracold atoms. Suppose an ensemble of atoms is 
prepared in the vibrational ground state of a harmonic trap. 
A Raman pulse with space-dependent Rabi couplings is ap- 
plied, as suggested in ifTTl . to put the atom in a superposition 
of dark states, after which the trap is switched off to allow bal- 
listic expansion. The initial wavepacket can be approximated 
by a Gaussian function g(k\G) = where d 

is the oscillator length of the trap, and k^^^ is a momentum 
boost (which is zero for the case of a stationary trap). For this 
wavepacket, the expectation values of y and z oscillate as 
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Due to the induced Born-Huang field, ZB will occur in this 
system (unlike its condensed matter counterparts |7l [H O) 
even if the wavepacket has an initial zero group velocity. Sup- 
posing ^^Rb atoms, we take the lattice wavenumber to be 
^0 = (27r/820) nm~^ and a Gaussian with = and width 
kod = 16.2, corresponding to the ground state of a trap with 
trap frequency 112 Hz ifTSl . Fig. |3] shows that a pronounced 
oscillation would occur in the z direction before damping out. 
A typical time scale of the ZB here would be /lis rather than 
the fs predicted in e.g. graphene and related systems 1 8j . 
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where we have defined L 
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shows that the ZB vanishes in the Abelian 
cases, 8 = or 8^ 1. 

It is useful to consider the limit d ^ oo^ where g(k',0) 
5(X — X^')), i.e., the initial wavepacket approaches a plane 
wave. The integrals in Eqs. ^TT\ become trivial, and we find 
that 
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In the opposite limit, J ^ 0, the ZB vanishes, and for inter- 
mediate values the energy spread causes damping, as can be 
shown analytically for bilayer graphene L8J . 



Figure 3: ZB for an atom with zero momentum spread (dashed), and 
for a momentum spread corresponding to the velocity spread of a 
cloud initially in a trap of frequency 112 Hz (solid). ZB oscillations 
for finite momentum spread damp out over time, but persist over sev- 
eral periods. 

We have shown that the mean position of the atom os- 
cillates sinusoidally. However, ZB can also be viewed in 
terms of state-resolved spatial dynamics. For the Gaussian 
initial wavepacket, it is not difficult to show that as the cen- 
ter of mass of the cloud is oscillating in the z direction, in 
the y direction, the wavepacket separates by spin, such that 
(ji,2(0) = ^hkot/m (see Fig. [4]). This spin separation, which 
coexists with the ZB, is a manifestation of the atomic spin 
Hall effect proposed in a different setup |[T6ll ; a related effect 
occurs in velocity- selective coherent population trapping 1 17 1. 

Fig. [5] shows the dynamics of the "spin-up" component 
of the wavepacket in a representative non- Abelian case. In 
essence, the effective magnetic field deflects spin-up and spin- 
down in opposite directions. As the two wavepackets sepa- 
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Figure 4: (Color online) Spin separation and ZB for 8 = 0(0.2)1, as 
indicated on the right colorbar. For each value of 8 the dashed tra- 
jectory corresponds to the atom's mean position in "spin-up," while 
the solid trajectory corresponds to the mean position of the atom in 
"spin-down;" open and closed circles indicate the respective ends of 
these trajectories. In the Abelian cases, 8 = and 1, the trajectories 
are straight lines; the trajectory for 8 = is the vertical line y = 0. 



where P is the center-of-mass momentum and £ = ^^^^ is 
the total spin. The Hamiltonian for the center of mass is thus 
again of a single-particle Dirac form (with a higher spin), and 
in a dilute cloud of N particles with two-body interactions, the 
center of mass of the cloud undergoes ZB. 

In this paper, we have examined the dynamics of an atom 
in a tripod level scheme on an optical lattice; this common 
experimental setup gives rise to a non- Abelian gauge field 
which is isomorphic to the spin-orbit interaction in 2D elec- 
tron gases. The idea of "atomtronics," or engineering atomic 
versions of semiconductor devices, has generated recent in- 
terest (TSj [191. The prospect of engineering artificial spin- 
orbit couplings suggests the possibility of atomic "spintron- 
ics," for example engineering atomic counterparts of devices 
such as the Datta-Das transistor 1201 , which have yet to be 
successfully realized with electrons. The tripod system here 
exhibits atomic ZB, with an amplitude many orders of magni- 
tude larger than that offered by Dirac electrons or recently dis- 
cussed condensed matter systems. We believe it is a promising 
candidate for the experimental observation of ZB. 



rate, the coupling between the components results in oscil- 
lating "tails" on the wavepackets, giving rise to ZB. The ZB 
decays as the wavepackets separate and cease to interfere. 




Figure 5: As seen in Fig. |4] the wavepacket separates by internal state 
as it jitters. This figure shows the time-evolved probability distribu- 
tion of the "spin-up" component of the wavefunction for an initial 
Gaussian wavepacket (8 = 8/?, kod = 5, hk^t /m=\0). The black dot 
indicates the mean position, which has drifted to the right. The mix- 
ing between internal states gives rise to "tails" on the wavepacket, 
resulting in ZB. The ZB damps as the internal states separate. 

We note that the Hamiltonian for N particles in the non- 
Abelian gauge field, interacting via two-body interactions, 
separates in center-of-mass and relative coordinates, R and p 
respectively. The Hamiltonian is then H = HcM^Hp, where 
is a function only of the relative coordinates p, and 
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